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Abstract. 

We calibrate the distorted wave Born approximation (DWBA) for electron impact 
excitation processes empirically. Differential cross sections (DCS) for the excitation of 
the 2p 5 3s, 2p 5 3p, 2p 5 As, and 2p 5 4p configurations of Ne and the 3p 5 4s and 3p 5 4p 
configurations of Ar by electron impact are calculated using DWBA for incident 
energies between 20 and 100 eV. The calculated results are compared with the absolute 
experimental measurements and other theoretical results. We found that the structure 
of the DCS can be well reproduced by the DWBA model while the magnitude is 
overestimated for most cases considered here. The differences in magnitude between 
DWBA and experiment are used to test the calibration of DWBA such that the DWBA 
can be used to describe laser-induced electron impact excitation processes. These 
processes are involved in the non-sequential double ionization of atoms in strong laser 
fields. 
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1. Introduction 

The process of electron impact excitation of atoms and ions is one of the most basic and 
important processes in atomic physics. Numerous theoretical methods have been used 
for the process calculations, including distorted wave Born approximation (DWBA) p] , 
second-order distorted wave model [2], i?-matrix method [3], and convergent close- 
coupling (CCC) calculations [3], among which the DWBA is the simplest. The 
sophisticated theoretical models, such as CCC and i?-matrix method, are capable of 
reproducing accurate angular differential cross sections (DCS), as well as the absolute 
magnitude. They are more suitable for low incident energies. For higher energies, both 
the integrated and differential cross sections predicted by DWBA are fairly accurate. 
However, it has been well recognized that, at low energies, the total cross sections 
(TCS) predicted by the DWBA signicantly exceed the experimental values. Ideally, one 
would use the i?-matrix approach for low energies, the DWBA for high energies, but at 
intermediate energies, neither method is efficient if a large amount of data is needed. 

The purpose of this work is to correct (by renormalization) the DWBA predictions 
empirically such that the DCS calculated from DWBA can be used for low collision 
energies. We will use the empirical formula proposed by Tong et al. [5] for the total 
excitation cross sections. Our ultimate objective is to apply the calibrated DWBA (C- 
DWBA) to simulate the correlated momentum distributions in nonsequential double 
ionization (NSDI) of atoms in strong laser fields. 

The process of NSDI of atoms in linearly polarized laser pulses is one of the most 
interesting and challenging topics in strong field physics. In NSDI, one electron that is 
first released near the maximum of the oscillating electric field may be driven back to 
revisit the parent ion when the electric field is near zero. When the returning electron 
collides with the parent ion with energies above the ionization threshold, it may kick 
out another bound electron, resulting in an (e, 2e)-like process. The returning electron 
may also excite the bound electron to a higher excited state which is subsequently 
tunnel ionized when the electric field increases again. Since the year of 2000, complete 
experimental measurements on the full momentum vectors of the two outgoing electrons 
along the direction of polarization of the laser pulse have become available [7J, and a 
number of theoretical studies have also been carried out. 

Recently, Chen et al. [8] have developed a quantitative rescattering (QRS) theory 
which has been applied to various rescattering processes induced by short intense laser 
pulses [9j [TOj [TTJ [121 H3J EE3] • The significant advantage of the QRS theory is that it treats 
the rescattering processes in the laser field as laser-/ree scattering processes, where the 
laser-induced returning electrons are described by a wavepacket. The QRS enables us 
to simulate two-dimensional correlated momentum distributions for NSDI quantitatively 
by calculating the triple differential cross sections (TDCS) for (e, 2e) [13] and the DCS 
for electron impact excitation of ions [H]. However, to obtain the correlated momentum 
spectra that can be compared with experimental measurements, one needs to evaluate 
the TDCS for (e, 2e) and the DCS for excitation for all possible momenta of the returning 
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electrons. For NSDI of atoms in strong laser pulses, the highest energy, Ef 13 -*, of the 
returning (incident) electron is determined by the laser field, which is less than 100 eV 
for typical 800 nm lasers. To simulate the correlated momentum distributions for NSDI, 
the DCS's for electron impact excitation of the parent ion at all incident energies from 
threshold to £'? nax are needed. To reduce the computational effort a simple and efficient 
theoretical model is desirable. In this work, we develop the C-DWBA for this purpose. 

The organization of this paper is as follows: In section II, the basic theory of 
DWBA for electron impact excitation is presented and the method to calibrate DWBA 
is introduced. In section III, the DCS of DWBA for electron impact excitation of Ne and 
Ar at incident energies below 100 eV are normalized and compared with the absolute 
experimental data. The normalization factors are then used to test the calibration for 
DWBA. 

Atomic units are used in this paper unless otherwise specified. 
2. Theory 

In this section, we present the DWBA theory on electron impact excitation of atoms and 
the method to calibrate the DWBA at low energies. The formulas presented here are 
generic and therefore can be easily applied to the processes of electron impact excitation 
of ions which are involved in NSDI. 

2.1. DWBA 

Suppose we have an electron with momentum kj which collides with an atom A, after 
the collision, the scattered electron has momentum k/, and one bound electron in atom 
A is excited to a higher energy bound state. In the frozen core approximation, the 
"exact" Hamiltonian for the whole system is 

H = + V A+ (n) - + V A+ (r 2 ) + — . (1) 

Z Z r 12 

where ri and r 2 are the position vectors for the projectile and the bound state 

electron with respect to the nucleus, respectively. This Hamiltonian can be rewritten 

approximately as 

Hj = + Ujin) - l -V\ + V A+ (r 2 ) (j = i, f). (2) 

In this equation, Ui (Uf) is the distorting potential used to calculate the initial 
(final) state wave function x k . (xu. f ) for the projectile. In the distorted wave Born 
approximation, the direct transition amplitude for excitation from an initial state to 
a final state \&/ is expressed by 

/=(xir / (l)« r /(2)|V i |* i (2)x£(l)>, (3) 
where V{ is the perturbation interaction, 

v l = H-H l = — + v A+ ( n ) - Uiin). (4) 
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In Eq. ([3]), the initial and final state wave functions for the projectile satisfy the 
differential equation 

and the bound state wave functions are eigenfunctions of the equation 
1. 



Xk,(n) = o (j = ;,/), (5) 



*i(r 3 ) = (j = *,/), (6) 



where ej (j = i, f) are the corresponding eigenenergies of the initial and final bound 
states which can be expressed as 

* j (r)=ip NjLj (r)Y L . M .(r) (j = ij). (7) 

The exchange scattering amplitude is given by 

g = ^ f (i)x^m\H2)xi(i))- (8) 

Finally, the differential cross section for electron impact excitation is given by 
da 4 kf 1 



N(27rY 



dtt hi 2Li + 1 

x £ E (\\f - 9\ 2 + \\f + 9\ 2 ) ■ 0) 

Mi=-LiM f =-L f v ^ ^ 7 

The prefactor N in Eq. ([9]) denotes the number of electrons in the subshell from which 
one electron is excited. 

The distorting potentials, Ui and Uf, used in Eq. fl5]) to calculate the wave functions 
for the projectile in the initial and final states, respectively, are not determined directly 
by the formalism. Here, we use static potentials which take the form as 

Ufa) = V A+ {n) + / rfr 2 ^M£ (j = ij). (10) 

As shown previously, V^+(r) in Eq. ([TO]) is the atomic potential used to evaluate 
eigenstate wave functions of the bound state electron. Here we use the effective potential 
from Tong and Lin [15] based on single active electron approximation, which is given by 

V A +{r) = , 11 

r 

where the parameters a^, as given explicitly in table 1 in Tong and Lin [15], are obtained 
by fitting the calculated binding energies from this potential to the experimental ones 
of the ground state and the first few excited states of the target atom. 

2.2. Calibration of DWBA 

The overestimate of DWBA on DCS can be corrected by using the empirical method 
proposed by Tong et al. [5] to evaluate the total cross sections for electron impact 
excitation: 

= <*^* {AE - )IE 'f (||) • (12) 
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where 



f(x) = - plnx +5 



\nx 



(13) 



x \ xj 



X 



In Eq. fll2p . AE is the excitation energy for a given transition, e is the eigenenergy of 
the corresponding excited state. The original formula given by Tong et al. does not 
have the prefactor a which is added in the present work to ensure that it reproduces the 
same cross sections as those from DWBA at high energies. The parameters in Eq. ({TBI) 
have been obtained initially by fitting to the convergent-close coupling excitation cross 
sections for hydrogen and He + . Explicitly, these parameters are (3 = 0.7638, 7 = 1.1759, 
and 5 = 0.6706. 

The total cross section of DWBA at fixed incident energy E t = kf/2 can be obtained 
from Eq. (Q by 



By matching the total cross sections from Eq. (fT2l) with those from Eq. (1141) at high 
incident energies, say Ei = 500 eV, one obtains the prefactor a in Eq. ffl~2]) for excitation 
of each configuration. To calibrate the DWBA at low energies, we define a scaling factor 



It is the the scaling factor in Eq. ( IT5|) that should be used to normalize the differential 
cross sections of DWBA at different incident energies. 

3. Results and discussions 

The perturbative nature of DWBA makes it overestimate electron impact excitation 
cross sections of ions at low energies. To calibrate the DWBA theory, one should 
compare its predictions to accurate theoretical results or absolute experimental 
measurements. Unfortunately, neither are easily available for atomic and molecular 
ions. Thus we use neutral Ne and Ar atoms for the calibration. 

In Figures 1-4, the DCS's for the excitation of the 2p 5 3s, 2p 5 3p, 2p 5 As and 2p 5 Ap 
configurations of Ne for incident energies below 100 eV from DWBA are compared with 
the absolute experimental data [161 Ej- For the 2p 5 3s configuration, the DCS's from 
the i?-matrix theory are also plotted for incident energies of 30 and 25 eV. It can be 
seen from Figs. 1-4 that the DWBA overestimates the DCS's for all the cases considered 
here. To get best overall agreement, different normalization factors are assigned to the 
DWBA for different configurations and different incident energies. For incident energies 
below 30 eV, the DCS's of DWBA could be 2-7 times higher than the experimental 
measurements. To see the distorting effect from the DWBA, the results of plane wave 
Born approximation (PWBA), in which plane waves are used to describe the projectile 
electron in both the initial and final states, are also displayed for comparison. For 
scattering angles greater than 30°, one can see that the PWBA fails completely in 
predicting the angular distributions even for incident energy of 100 eV. In contrast, the 




(14) 



C(Ei) = a-T ong (Ei)/aBWBA(Ei). 



(15) 
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Figure 1. DCS for the excitation of the 2p 5 3s configuration of Ne by electron impact 
at incident energies of (a) 100 eV, (b) 50 cV, (c) 30 eV, and (d) 25 eV. The absolute 
experimental measurements are from Register et al. |16j and Khakoo et al. |17j . For 
incident energies of 30 and 25 eV, the results of i?-Matrix from Khakoo et al. [T7J are 
also plotted for comparison. 
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Figure 2. Same as Fig. 1 but for configuration of 2p 5 3p at incident energies of (a) 50 
eV and (b) 30 eV. 

enhanced DCS's for backward scattering observed in experiment are well reproduced by 
the DWBA. 

In Figs. 5 and 6, we show similar comparison for the excitation of 3p 5 4s and 3p 5 4p 
configurations of Ar. The experimental measurements were performed by Chutjian and 
Cartwright [18]. Compared to the excitation of Ne, the DCS of Ar have more structures. 
For example, for 3p 5 4s at 100 eV and 50 eV, as shown in Figs. 5(a,b), in addition to the 
rapid slope change around 25°, extra minima were observed in experiment which are 
reproduced by the DWBA. For 3p 5 4p at 100 eV and 50 eV, as shown in Figs. 6(a,b), 
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Figure 3. Same as Fig. 1 but for the 2p 5 4s configuration at incident energies of (a) 
50 eV and (b) 30 eV. 
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Figure 4. Same as Fig. 1 but for the 2p 5 Ap configuration at incident energies of (a) 
50 eV and (b) 30 eV. 
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Figure 5. DCS for the excitation of the 3p 5 4s configuration of Ar by electron impact 
at incident energies of (a) 100 eV, (b) 50 eV, (c) 30 eV and (d) 20 eV. The absolute 
experimental measurements are from Chutjian and Cartwright |18j . 
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Figure 6. Same as Fig. 5 but for the 3p 5 4p configuration at incident energies of (a) 
100 eV, (b) 50 eV, (c) 30 eV, and (d) 20 cV. 



DWBA predicts triple minima in the DCS. They were observed in experiment as well 
despite that the backscattering is overestimated by the DWBA. This might indicate 
that the distorting potential used in the calculations needs to be improved. For lower 
incident energies of 30 eV and 20 eV, the agreement between DWBA and experiment 
for 3p 5 Ap of Ar can not be regarded as satisfactory. However, the main feature can still 
be predicted by the DWBA. It should be noted that, for the excitation of 3p 5 4p of Ar 
at 20 eV (see Fig. 6(d)), the yield from DWBA exceeds the experimental value by an 
even greater margin than the PWBA result. For this case, the DWBA predicts the DCS 
which is about 40 times higher than experiment. 

The empirical formula, Eq. (I12p . has already been used to calculate the total 
ionization yield of Ar in NSDI as a function of the peak intensity for a linearly polarized 
laser pulse by Micheau et al. [12] . 

To obtain the scaling factor C(Ei) for electron-Ne excitation, we calculate the TCS's 
using DWBA and those using Eq. ( 1121) for all the four configurations considered here 
at incident energies from threshold up to 500 eV. These TCS's for incident energies 
below 120 eV are plotted in Fig. 7 referring to the left vertical axis. It can be seen 
that the difference in the magnitude of TCS between DWBA and Tong et al increases 
with decreasing incident energy. This difference is indicated by the the scaling factor 
C(Ei) which is also plotted in Fig. 7, referring to the right vertical axis. To see the 
accuracy of the scaling factor, the normalization factors used in Figs. 1-4 for DWBA 
to obtain the best overall agreement with the experimental DCS's are displayed for 
comparison. One can see that all the normalization factors used in Figs. 1-4 agree well 
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Figure 7. Total cross sections (left vertical axis) and normalization factors of DWBA 
(right vertical axis) for electron impact excitation of Ne from 2p 6 to (a) 2p 5 3s, (b) 
2p 5 3p, (c) 2p 5 4s, and (d) 2p 5 4p. Solid curve, total cross sections of DWBA; Dotted 
curve, total cross sections calculated using the empirical formula of Tong et al. [5]; 
Chain curve, scaling factor C(Ei); Solid circles, normalization factors used in Figs. 1-4 
for DWBA to obtain the best overall agreement with experiment. 
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Figure 8. Same as Fig. 7 but for Ar from 3p 6 to (a) 3p 5 4s and (b) 3p 5 4p. The solid 
circles are the normalization factors used in Figs. 5 and 6 for DWBA to obtain the 
best overall agreement with experiment. 



with those predicted by the scaling factor C(Ei). In Fig. 8, similar comparisons for Ar 
are shown. The good agreement of C(Ei) with the normalization factors used in Figs. 5 
and 6 confirms again the validity of the calibration method. 

In conclusion, we proposed a method to calibrate the DCS from DWBA for electron 
impact excitation of atoms at low energies. The method will be applied to simulate the 
correlated electron momentum spectra for NSDI of atoms in a strong laser field, in 
which electron impact excitation of the ions is involved. This work paves the way for 
theoretical study, based on the QRS model, on NSDI of atoms in strong laser pulse. 
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